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ABSTRACT 


This study considers the effects of injection of 
a nonreacting gas (taken to be nitrogen) on the heat 
and mass transfer of a porous, chemically reacting 
surface (graphite) at hypersonic speeds. Heat and 
mass transfer relations are obtained for the carbon 
surface. 

Numerical results are presented for the case 
of the sphere stagnation point at hypersonic speeds. 
Reductions in carbon surface mass loss and aerody- 
namic heating as a result of nitrogen injection are 
demonstrated. Comparisons with more rigorous 
digital computer solutions for zero nitrogen injection 
show the same qualitative trends and numerical dif- 
ferences of about 10 percent. 
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AN APPROXIMATE ANALYSIS OF THE EFFECTS OF NITROGEN 


INJECTION INTO THE BOUNDARY LAYER ON A GRAPHITE 

SURFACE IN A HIGH- TEMPERATURE DISSOCIATED AIRFLOW 

By Kenneth C. Weston 
Manned Spacecraft Center 

SUMMARY 


This study considers the effects of injection of a nonreacting gas (taken to be 
nitrogen) on the heat and mass transfer of a porous, chemically reacting surface 
(graphite) at hypersonic speeds. All chemical reactions are considered to occur in a 
reacting sheet adjacent to the wall. Equations for the determination of species mass 
fractions at the wall are derived under the assumptions of a frozen boundary layer and 
chemical equilibrium in the flame sheet. Schmidt and Prandtl numbers are taken as 
unity. Heat and mass transfer relations are obtained for the carbon surface. 

Numerical results are presented for the case of the sphere stagnation point at 
hypersonic speeds. Reductions in carbon surface mass loss and aerodynamic heating 
as a result of nitrogen injection are demonstrated. Comparisons with more rigorous 
digital computer solutions for zero nitrogen injection show the same qualitative trends 
and numerical differences of about 10 percent. 


INTRODUCTION 


Graphite is considered frequently as a heat -protection material for the leading 
edges and other high heating-rate areas of lifting reentry vehicles. Its ability to sus- 
tain high temperatures and its strength at elevated temperatures make it particularly 
suitable for such applications. These unusual qualities have motivated considerable 
analytical study of the ablative behavior of graphite. Moreover, its elemental nature 
offers substantial analytical simplifications relative to chemically complex organic 
ablators. The present work utilizes this simplicity to examine the effects on heat and 
mass transfer of injection of nonreacting nitrogen through porous graphite. 

A number of investigators have analyzed the behavior of the flow over an ablating 
substance within the framework of boundary -layer theory. In 1957, Dennison and 
Dooley (ref. 1) treated the laminar boundary layer with combustion using a model very 
much like that employed in the present study. This model assumes all chemical reac- 
tions to occur at a surface to which oxygen and fuel are transported by convection and 
diffusion. Chemical equilibrium is assumed to exist at this surface; combustion 
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products diffuse in both directions from the surface. As a further simplification, ref- 
erence 1 utilized constant property boundary -layer solutions for Lewis and Prandtl 
numbers of unity. 

In 1962 and 1964, S. Scala (refs. 2 and 3) presented the results of sophisticated 
calculations of the ablation of graphite. These digital computer solutions of the re- 
acting gas boundary layer included consideration of nine chemical components and six 
chemical reactions in the boundary layer. The calculations involved simultaneous 
integration of coupled boundary -layer equations including calculations of transport 
properties based on multicomponent diffusion rather than the frequently employed bi- 
nary model. These calculations show distributed reaction zones within the boundary 
layer, and are the most detailed results presently available for graphite ablation. 

In reference 4, Lees concluded that, for a Lewis number of unity, reaction rates 
in the boundary layer play a secondary role in determining the mass fractions of chem- 
ical species at the surface and transport properties. His analysis also suggests that 
combustion may be beneficial in reducing heat transfer if the heat of sublimation is 
comparable to the heat reaction. Lees further stated that these conclusions are inde- 
pendent of the model of chemical reaction assumed. 

Dorrance outlined a method in reference 5 for the analysis of the reacting lami- 
nar boundary layer on a carbon surface. His analysis considered all reactions to 
occur at the wall, with the boundary layer treated as nonreacting. The present work 
extends this analysis, for the case of Prandtl and Schmidt numbers of unity, to the 
case of nitrogen injected into a porous carbon surface. 

This work is based on research documented in a Master’s Thesis at Rice Univer- 
sity in 1965. Discussions with Dr. Frederic A. Wierum of Rice University are grate- 
fully acknowledged. 


SYMBOLS 


B mass transfer parameter 

C species mass fractions 

C elemental mass fractions 

C„ Stanton number 

H 

0 ^ specific heat at constant pressure 

D diffusion coefficient 

f modified stream function (see equations (59), (60), and (61)) 

G function defined by equation (74) 
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g dimensionless enthalpy, y- 

e 

(g) gaseous phase 

h local static enthalpy 

h ,, effective heat of ablation 

ert 

h. enthalpy of the i'th specie 

In 0 heat of reaction of the i'th specie 

I stagnation enthalpy 

Jj, Jg, Jg quantities defined in equations (38) to (40) 

Kg,K^ equilibrium constants defined in equations (1) to (4) 


k 

L 

L 

\ 

a 

M 

M 


M* 


M * 
e 


thermal conductivity 
pDc 

Lewis number, — s- 


P 


heat of sublimation of graphite 


density -viscosity ratio defined in equation (63) 
molecular weight 

mean molecular weight defined in equation (81) 
Mach number 

molecular weight of undissociated air 


Mach number at edge of boundary layer 


rii£ mass flux of carbon 

ii> total mass flux from equilibrium layer defined in equation (10) 

E 


mass flux of nitrogen 


Prandtl number, 


r P 

IT 

pressure 
heat flux 

energy absorbed in steady ablation 
radiative heat flux from surface 
heat flux to graphite surface 

aerodynamic heat flux with mass injection effects included 

sphere radius of curvature 
universal gas constant 
gas constant for i'th specie 

radius normal to body axis of symmetry extending to body surface 
Schmidt number, p/pD 
coordinate along body surface 
solid phase 

boundary -layer coordinate defined by equation (58) 
temperature 

reference temperature defined by equation (79) 
time 

local velocity parallel to wall 
local velocity normal to wall 
rate of production of the i'th specie 



-1 


, equation (30) 



x i 


y 


z i 


6 

TJ 

T/*,7J** 

X 

P 

V 

P 

T 

* 


mole fraction of the i'th specie 
coordinate normal to wall 

dimensionless mass fraction of i'th specie, ^ 

C i,e 

thickness of equilibrium layer 

boundary -layer coordinate defined in equation (57) 

dummy variables 

generalizing variable representing f, g, or 
viscosity coefficient 

index, v = 0 for planar flow, v = 1 for axisymmetric flow 
density 

thickness of graphite layer 

generalizing variable representing 1, P, or S 
stream function 


Subscripts: 

A air 

E condition in equilibrium layer 

e condition at outer edge of boundary layer 

f frozen 

i of i'th specie 

in initial 

M arithmetic mean value 

N. B. no mass injection 

stagnation point 


sp 



w 


conditions at wall 


OO 


ahead of bow shock 


ANALYSIS 


Analytical Model 

The ablation of graphite at hypersonic speeds occurs through sublimation and 
chemical reactions with the surrounding air. At wall temperatures below about 
1000° K the rate of mass loss is governed by the rate of oxidation of the graphite sur- 
face. At slightly higher temperatures a transition region precedes a third region in 
which the rate of surface recession is governed by the rate of diffusion of oxygen into 
the boundary layer. This "diffusion -controlled regime" prevails over a range of wall 
temperatures from about 1500° to about 3000° K at atmospheric pressure. The pres- 
ent analysis considers the effects of nitrogen injection on graphite ablation in this re- 
gime. At temperatures above 3000° K the rate of sublimation of carbon exceeds the 
rate of diffusion of oxygen, and uncombined gaseous carbon molecules exist at the wall 
in significant amounts. In this regime nitrogen reactions with graphite yield sizable 
mass fractions of cyanogen. 

The present study considers the ablation of graphite for the diffusion -controlled 
regime in terms of a model in which the rate of diffusion of oxygen to the chemically 
active graphite surface is governed by mathematical solutions for a frozen boundary 
layer. The boundary layer is frequently termed frozen when the rate of generation of 
species in it is negligible. Such solutions have been successful in predicting heat 
transfer in dissociated-air boundary layers in which the wall catalyzes atom recombi- 
nation and the wall temperatures are well below the boundary -layer -edge temperature 
(refs. 4 to 6). It is therefore assumed that all chemical reactions occur in an infini- 
tesimal layer adjacent to the wall in which chemical equilibrium exists. This "equilib- 
rium layer" is a region where nitrogen and carbon are received from the ablating 
surface, chemical reactions occur with oxygen from the boundary layer, and the equi- 
librium mixture is supplied to the boundary layer. In view of the approximation im- 
posed by some of the necessary assumptions, it is reasonable to take advantage of the 
great simplifications made possible by assuming Prandtl, Schmidt, and, therefore, 
Lewis numbers to be unity. These assumptions are perhaps best justified by the re- 
markable agreement of the present calculations with the results of reference 3 which 
were obtained with state-of-the-art transport properties. 

The present analysis is limited to hypersonic conditions for which the species at 
the edge of the boundary layer are atomic or molecular. Moreover, the boundary - 
layer-edge composition is assumed to be made up of nitrogen and oxygen only. 
Boundary -layer -edge conditions are obtained from the normal shock conservation equa- 
tions, the assumption of isentropic and isoenthalpic flow along the body streamline, and 
the real air properties of reference 7. Absence of consideration of ions and electrons 
at the boundary -layer edge limits the applicability of this analysis to boundary -layer - 
edge temperatures below 10 000° K. This is well beyond the maximum equilibrium- 
shock-layer temperatures encountered, for example, during the reentry of current 
manned orbital vehicles. 
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The thermal decomposition of a graphite surface has been studied by many inves- 
tigators. While carbon molecules (atomic carbon C and higher mass carbon radicals 
such as C 0 ,CL, C., . . . C ) may sublime from a hot carbon surface, experimental 

evidence indicates that C g and C are most abundant over a wide range of tempera- 
tures (refs. 3 and 8). In the present analysis it is assumed that the sole product of 
carbon mass loss is gaseous C^. Dissociation of Cg to C and subsequent reactions 

with oxygen to produce carbon monoxide are allowed in the equilibrium layer. Forma- 
tion of carbon dioxide is not considered since it has been shown in reference 3 to be 
present only as a trace specie in the temperature range of interest. The analysis is 
therefore based on the assumed heterogeneous surface reaction 


C 3 (s) 



c 3 (g) 


and the following homogeneous reactions in the equilibrium layer 


C 3 (g) ^±3C(g) 
2C(g) +0 2 (g)^=^2CO(g) 


C(g) + 0(g)^=^C0(g) 

where (s) and (g) indicate solid and gaseous phases, respectively. 

The results of reference 3 show that nitrogen reactions are insignificant in the 
temperature range -of -interest. In this case nitrogen acts only as a diluent. More- 
over, nitrogen is present in the equilibrium layer in diatomic form only; the range of 
wall temperatures considered is too low for significant dissociation. 


Composition of the Equilibrium Layer 

The partial pressures of the species in the four chemical reactions listed may be 
related through the following expressions for the equilibrium constants: 


K.(T) = p c (1) 

1 °3 

K 2 ( T ) - p c Vp C3 (2) 
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(3) 


K 3< T > - PC0 2 A/C 2 

K 4 (T) = P C( yVo <4 > 

The partial pressure of Cg is given by Thorne and Winslow (ref. 8) in the vi- 
cinity of 2400° K as 


, „ „„„ 40 296 

log p c ^ = 9. 811 - ■ vf - 


(5) 


where T and p are specified in °K and atmospheres. Thus the partial pressure of 
Cg at the wall is determined when the wall temperature is specified. The partial pres- 
sures of the other species, except diatomic nitrogen Ng, may be expressed in terms 
of the partial pressure of carbon monoxide and the equilibrium constants as follows: 


>0 = (v\ 


K */ 3 

K 1 K 2 

(6) 

r -2/3 -2/3\ 2 

"1 K 2 j p CO 

(7) 

-1/3 -l/3\ 

1 K 2 j p CO 

(8) 


The pressure at the wall is related to the partial pressures of the above species 
and nitrogen by Dalton’s law of partial pressure 


p N 2 + P 0 2 + p O + P C + p Cg + P CO P E P e P w 


(9) 


The latter equality follows from the usual assumption of constancy of pressure 
across a boundary layer. To obtain the expression for the partial pressure of carbon 
monoxide necessary to determine partial pressures of the other species in the equi- 
librium layer, the convection and diffusion of species and mass injection at the wall 
must be considered. 
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Mass Balance 


Consider the flux of mass through the equilibrium layer at the wall as shown in 
the following sketch. In general, nitrogen and carbon are uniformly injected into the 


e 

E 

w 


m 


E 


m 


N, w 


m 


C, w 


Boundary layer 


Equilibrium layer 


equilibrium layer and a net mass flux m E from the equilibrium layer into the bound- 
ary layer is defined. Then for steady flow 


m. 


= m„ + m_, 

N, w C, w 


(10) 


Expressions may also be written for the accumulation of mass of each of the ele- 
ments C, N, and monatomic oxygen O in the infinitesimal equilibrium layer. For a 
steady flow of each of the elemental mass fractions, the net rate of storage of each 
element must be zero. The following relations must then be satisfied: 


dC C _ f 0C C 

6p -df = *C, w - Vc,E + \^~sr) E = 0 


( 11 ) 


dC N _ l 3C n 

6p lir = a n,w- A E C N,E + l pD -aT J =0 


( 12 ) 


dC o _ / 3C o) 

6p ~ST “ " m E C 0, E + ^ pD 1yy E = 0 


( 13 ) 
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For later reference, equations (11) to (13) may be written 


_ 5C C 

m_, + m XT = m. T + - I pD -=■ — 

C, w N, w N, w E C, E \ dy 


E 


ac 


“W = “A th ' PD 


N 


N,w E N, E V " dy 


dc ( 

0 = m E C 0, E - ( PD 1F 


IE 


(14) 


(15) 


(16) 


It is shown subsequently that when the Prandtl, Schmidt, and Lewis numbers are 
unity, the velocity, enthalpy, and elemental mass fraction profiles are linearly related. 
Therefore, 


where 


Then 


c = c E 





dT j + ^ c i h i° + u V 2 


ac _ c e " c e ai 
ay " i e -i E ay 


(17) 


and the diffusion term in equations (14) to (16) may be expressed in terms of the Stanton 
number 



' L if 
C P 




P u C TT 
e e H 


(18) 
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where L = 1 and the Stanton number C„ is given by 

n 


/k af 

I c dy 

\P > 


l 


p u C„ = T - 
e e H I - I T7 , 
e E 


The diffusion terms in equations (14) to (16) may be eliminated by using equa- 
tion (18). The resulting expressions are 


n E " m N, w + m E C C, E ‘ ( C C, e ' C C, e) P e U e C H 
“N, w " *E C N, E - (S, e ' U N, E ) P e U e C H 
0 = m E C 0, E " ( C 0, e " C 0, e) P e U e C H 


(19) 

( 20 ) 
( 21 ) 


where equation (10) has been used to eliminate the injection rate of carbon at the wall. 
A mass injection parameter B is defined as 


m 


B = 


E 


p e u e C H 


( 22 ) 


Equations (19) to (21) become 


m 


N, w 


<B + ,)C CiE = C Cie + B-^ 


(23) 


+ TP “ C-KT Q + 


m 


N, w 


'N, E N, e p u C„ 
’ e e H 


(24) 


(B + l)U 0>E = ^o,e 


(25) 
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Throughout the following pages the subscript w is omitted from the wall 

ih, T 

N w 

nitrogen-injection parameter — b — with the understanding that nitrogen injection 

" P U v T r 

*e e H 

through the wall is indicated. 

Equations (23) to (25) represent the mass balances on the equilibrium layer in 
terms of the elemental mass fractions. These expressions may be rewritten in terms 
of species mass fractions. In so doing, the mass fractions of carbon at the edge of the 
boundary layer and the mass fraction of monatomic nitrogen at the wall are set equal 
to zero. 


(B + 1) 


C C„ E + C C 


'■Tcoj 


CO, E 


m 


= B - 


N 


P u C IT 
e e H 


(26) 


(B 


< B + » ( C N 2 , e) = ( C N 2 , e + C N, e) + pu^C R 
+ ( C 0, E + C 0 2 , E + C CO, e) = ( C G, e + C 0 2 , ej 


(27) 


(28) 


It is desirable to write the mass fractions in the equilibrium layer in terms of 
partial pressures of the species. The partial pressure and the mass fraction of the 
i'th specie are related by 


C. = p.M.X 
i *i i 

where 



(29) 


(30) 
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Equations (26) to (28) become 


(B + 1) fM c ^p c ^ E + M c P C, E + M C* 




® + 11 ( M n/n 2 , E ) X E - C N 2 , e + C N, e + P e u”c H 
(B + 1) e + M 0 ^p 0 ^ E + MqP^ e ^X e = C Q e + C Q 2j 


All of the basic relations are now available to solve for the partial pressures of 
the species in the equilibrium layer. From equation (32) the partial pressure of dia- 
tomic nitrogen Ng is given by 


N 2’ E 


C N 2 , e + C N, e + p e u e C R 
(B+DXeM 


Equation (33) may be used to eliminate (B + 1)X F from equation (34) yielding 


N 2’ E 


C N 2 , e 4 C N, e + PgUgCjj J (^cfo, E + M 0 2 P C> 2 , E + M oPcO, e) 

M_ t /C« + C„ \ 

N 2 \ e ° 2 ’ e ) 
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Using equations (7) and (8) the partial pressure of nitrogen may be expressed in terms 
of the partial pressure of carbon monoxide 


m. 


’n 2 ,E 


1 C N 2 , a + C N, a + pTTC H )( M 0 K 4‘ lK r 1/3K 2 ' 1/3 Pc0 > E ) 

M„ tC 


N 2 ( C °, e + C 0 2 , e) 


nu 


l C N 2> a + C N, a + l^bo 2 K 3'\' 2/ \' 2/ \o, e' + “0*00, e) 


M N 2 ( C °, e + C 0 2> e) 


(36) 


Combining equation (36) with equations (1) and (6) to (9) yields a quadratic 
equation for p CQ £ 


m N ) 

M N,ro,e + t O„ e L 


P CO, E 


-1 -2/3 -2/3 

+ M q Kg Kj K 2 P co 

6 


1; e 2 ] + (k 3 -‘k 1 - 2 / V 2 / >co,e 2 

( K 4 ' lK r 1/3 K 2 ' 1/3 )p CO , E + k/ 7 V /3 + K 1 + PCO, E - Pe = 0 


(37) 


Let 



N 

% a + C N, e + 

*\(Co,a + <V) 


(38) 


14 



II 

? 


J 2 = 1 + k 4 -v v V 1/3 


_ K -1 K -2/3 -2/3 

J 3 " K 3 K 1 K 2 


Expressed in terms of the J, equation (37) becomes 


2 J 1 J 2 M 0 + J 2 


K 1 1 / \ 1 /Z +K 1 -p e 


? CO, E + J 1 J 3 M 0 +J 3 P CO, E + +j 3 

2 2 


= 0 


(39) 

(40) 


(41) 


The complete solution for the partial pressure of carbon monoxide may be shown as 


p CO, E 


1/2 


l( J 1 J 2 M 0 
° 2 




- K, - K/V^\ 

J 1 J 3 M 0 2 + J 3 j 


J 1 J 2 M Q + J 2 
J 1 J 3 M 0 2 + J 3 


(42) 


Over a wide range of wall temperature and pressure, the squared term in the 
radical of equation (42) is large compared with the other term. Thus numerical diffi- 
culties are encountered in attempting to evaluate p^Q E using equation (42). Exami- 
nation of the orders of magnitudes of the terms in equation (41) reveals that while 
2 

P CO E * s or< * er one ’ th e °ther two terms are many orders of magnitude larger. 

In this case p^Q E can be accurately calculated by neglecting p^Q E ^ in equa- 
tion (41). The following simpler expression is obtained 


CO, E 


p e ~ K 1 ~ K 1 K 2 
J 1 J 2 M 0 + J 2 


(43) 


This expression may be accurately employed for the determination of 

when 


P CO, E 


J 1 J 2 M 0 + J 2 


J 1 J 3 M 0„ 


+ J„ 


» P 


CO, E 


(44) 
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Using equation (43) the partial pressures of the species other than CO may be 
expressed in terms of J as 


„ . k - K i - Ki 1, 

P °2’ E ' 3 \ J 1 J 2 M 0 + 


1/3 l/3> 

K 2 


/ A - K i - k i 1Ak ' 

P 0, E ( J 2 " ^ J i J 2 M o + J : 


1/3 v 1/3 


_ 1/3 1/3 

P C, E “ K 1 K 2 


P C 3 , E " K 1 


P N 2 , E J 1 ( M 0 J 2 P C0, E + M 0 2 J 3 P CO, E ) 


_ P e " P CO, E ' P 0 2 , E ' P 0, E " P C, E ' P C y E 


(45) 

(46) 

(47) 

(48) 


(49) 


When the wall temperature, nitrogen injection parameter, and boimdary- layer - 
edge conditions are specified, J can be calculated from equations (38) to (40); the 
equilibrium -layer partial pressures can be evaluated from equations (43) and (45) 
to (49), and B is determined from equation (32). To evaluate the heat and mass 
fluxes at the graphite surface the theory of the boundary layer must be considered. 


Boundary Layer 

The present analysis of the boundary layer follows the treatment of Dorrance in 
reference 5. The basic equations for the reacting gas boimdary layer are presented. 


16 



Continuity: 


a 

da 



(50) 


S- momentum: 


3u 

pu^ + pv 


3u 

3y 




(51) 


Y-momentum: 


P = 


p e = p E 


(52) 


Energy: 


ai ai a I u 

PU-5Z + PV — = — 


3s 3y 3y 


91 . „ 4 Aiaufl a |?1 jL n v*. aCi 
ay + ' pj 2 dy J " 3y |^L " 7 pD £ h i 3y 


(53) 


Species: 


pu 


ac. ac. , / ac.\ 

~a^ + ^ “ 3 y = W ( pD ” 9 y j + *i 


(54) 


State: 


Pj = P,R,T 


(55) 


Enthalpy: 


-Hi c p,i dT+h i°V' 12 / 2 = h + u 


2 / 2 


(56) 
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Equations (50) to (56) apply to the flow of a laminar boundary layer over a plane 
(v = 0) or axisymmetric {v = 1) surface. The gas is assumed to be an ideal gas and the 
diffusion process is assumed to be represented by binary diffusion according to Fick’s 
law. Equations (50) to (56) are a set of nonlinear partial differential equations which 
are not generally solvable. A very useful technique for the solution of such equations 
is to transform the equations into ordinary differential equations. Solutions of this type 
are termed similar solutions. The transformation of equations (50) to (56) is derived 
in detail in reference 5. The equations are transformed into an rj, s coordinate sys- 
tem in which 


pur 1 ' rY 
„ _ e e 
’■wF 71 


/. 


o Pe 


dy 


(57) 


and 




2v , 

pur u ds 
ee e 


(58) 


A stream function which is the product of separate functions of rj and s is 
introduced 


i//(s, rj) = N(s) f(v) 

The function f (rj) is defined such that 

u = u g f ' (rj) 


(59) 


(60) 


The modified stream function at the wall f(0) is also related to the mass flow 
through the wall as 


f(0) = 


-(pv) E (2s) 17 2 


p u u r 
e e n e 


(61) 
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The dependent variable in the momentum equation becomes f, and & is elimi- 
nated. The resulting equation is then considered to be a function of r) only. The 
resulting momentum equation becomes 


(£f")’ +ff” 


u e ds [p 



0 


where primes denote differentiation with respect to 7? and 


(62) 


i = 


PP 


Ve 


(63) 


It should be noted that similar solutions may be obtained for equation (62) only if 
the third term is a constant or a function of tj only. This is the case for a stagnation 
point where 


n — UU 

2s e , , 

- — = a constant 

u e ds 


or for a flat plate where 



= 0 


If dimensionless enthalpy and species functions are defined as 


g(v) = 



Z.(rj) = 



(64) 

(65) 


19 



transformations of the energy and species equations may be written 



and 


(1 z i’)' * iZ i * 


2§f ’ Z. dC i 
1 e_ 

C. ds 
1 , e 


2sw. 

l 


2 2Vr, 

pp u p r C. 
e e r e i, e 


(67) 


The boundary conditions on these equations are 


f(°) = f E 


-(pv) E (2s) 


1/2 


p eW 


f’(0) = 0 

f*(°o) = i 

C i E 

z.(o) = z. -Ji* 

1 X E c i, e 

Z.(^) = 1 

g(o) = g E 

g(“) = 1 


At this point it is desirable to attempt to simplify equations (62), (66), and (67). 
Lees (ref. 9) has pointed out that, in the presence of extreme cooling, the third term 
in equation (62) has only a weak effect on solutions of the stagnation-point problem. 
Therefore, it is assumed that this term is negligible compared with the first and sec- 


du 

ond terms. On a flat plate with no pressure gradient -g=- = 0 and the third term is 


identically zero. In the present case the momentum equation is adequately represented 
by 


Uf")' +ff" = 0 


( 68 ) 
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The external flow at the edge of the boundary layer is considered to be isoener- 

dl 

getic. This leads to = 0 and the removal of the first term on the right-hand side 

of equation (66). The external flow at the boundary- layer edge is evaluated from equi- 
librium calculations. The species mass fractions are assumed to be frozen at these 
equilibrium values in the case of both the flat plate and stagnation point. This leads to 
the elimination of the first term on the right-hand side of equation (67). As discussed 
in the analysis section, this report considers the case in which no chemical reactions 
occur in the boundary layer and all reactions are at equilibrium in an infinitesimal 
layer at the wall. The rate of generation of the various chemical species in the bound- 
ary layer is then zero; hence w. = 0. Finally, with the assumption of Prandtl, 

Schmidt, and Lewis numbers equal to unity, equations (66) and (67) may be written 


+ 

crq^ 

II 

o 

(69) 

(r z i ) + tz i - 0 

(70) 


) 


It is evident that equations (68) to (70) are of identical form and may be written 


+ fX' = 0 (71) 


P, or S. It is clear that f’, g, and Zj 

are linearly related when P = S = 1 since they all satisfy the same differential equa- 
tion. The functions f’ and g are also linearly proportional to the elemental mass 
fractions because the latter are proportional to Z^. The boundary conditions associ- 
ated with equation (71) are 


w 


where X is f’, g, or Z. and <P is 1, 


X(0) = X E , V(0) = X E ', X(«) - 1 


Equation (71) is readily integrated with (f> and £ as constants to yield 


X* = X' (0) exp 



(72) 
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and 


Let 


X = X(0) + X’(0) J ^ exp^- j £f dn**jdn* 


G(v, (t>) = fj j exp^- j f d ti**J dq* 


so that using equation (73) and the boundary conditions X' (0) and X(rj), 


x.(o)=-L_*m. 


0 G(°°, <p) 


and 


Mo) = MO) + [1 - MO)] 

It is evident from equation (75) that if G(°°, 0) is known, g' (0) and 
therefore, the convective heat transfer can be determined, since 


g’(0) = 


1 - g(o) 

| G(°°, P) 


and 


Z.’(0) 


1 - Z.(0) 
J-G(°°, S) 


(73) 

(74) 

(75) 

(76) 

;.'(0) and, 

(77a) 

(77b) 
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Figure 1. - Heat transfer parameter. 


Constant property solutions of equa- 
tion (71) are available in references 10 and 
11. Results from these references are pre- 
sented in figure 1 which relates G(°°, 4>) to 
f(0). The mass injection parameter B may 
be related to these functions using 


«°> * < 78 > 


which follows from equations (22), (61), 
and (92). 

Thus, given P = S = 1, £, and B, 
the heat transfer parameters g’ and Z' 
are readily obtained from figure 1 and equa- 
tion (77). 


Property Determination 

Dorrance (ref. 5) recommends the use 
of the reference temperature in evaluating 
6 . . The concept of a reference temperature 
has been utilized for some time in simula- 
ting the results of variable property solu- 
tions of the boundary- layer equations. 

Eckert (ref. 12) defines a reference tem- 
perature as follows 


= 0. 5T + 0.038TM* 
e e e 


+ 0. 5T 


E 


(79) 


The evaluation of i using the reference temperature is straightforward for air 
or pure gases whose properties are well known. In the present case, however, the 
question arises as to the proper method of evaluating £ for a mixture of gases. The 
following approach is adopted: 

(1) The viscosity is taken as the viscosity of air at the temperature T + using 
the Sutherland equation 




= 3.06 X 10 


t 3/2 
T + 112 


(80) 
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, _ . . . . lb-sec 

where T is in K and p in « — . 

ft 2 

(2) The molecular weight M is calculated using the arithmetic mean of the spe- 
cies mass fractions at the wall and the boundary-layer edge as 


M = 



where 


C. 

h 


M “ 


i, e 



(81) 


(82) 


(3) The density may then be calculated from the perfect gas law 



(83) 


This approximate approach for evaluating £ (t + ) was compared with a more de- 
tailed calculation for the stagnation-point boundary layer in which species mass frac- 
tions were obtained at the point in the frozen boundary layer where the temperature is 

T + . These species mass fractions were then used to evaluate the molecular weight 
and viscosity. No significant difference was obtained between the results of this tedi- 
ous calculation and the simple method presented herein. In view of existing uncertain- 
ties in parameters necessary for detailed calculations of transport properties at high 
temperatures and the approximate nature of the present analysis, it is thought that 
more sophisticated approaches are unwarranted. 


Convective Heat Transfer 

Heat transfer at the carbon surface is assumed to be the result of thermal con- 
duction and the transport of chemical energy which is released or absorbed at the wall 
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II 

f 


because of chemical reaction. The wall heat flux is assumed proportional to the gra- 
dients of temperature and species mass fraction at the edge of the equilibrium layer 


^w 


= K 


9T 

3y 


ac. 

+ P D Z>i T? 


E 


(84) 


The temperature gradient may be related to the enthalpy gradient by 

rT 


(4 - (? I i dT S -w , f * ZV ^ 


The frozen specific heat c^ ^ may now be introduced where 

c t = Y\ C.c 
p, f V i P> i 


E 


(85) 


( 86 ) 


so that 


where 


(a- 


'p, f 


„ ac.\ 
ai ya ii 

^3y 4" i ay j 


E 


(87) 


„ /-T ac. _ ac. / rT \ ac. ac. 

2/ c . dT^ + ^h .°-=-i = V[ / c . dT + h ° ) = Vh. (88) 

i P) 9y i 1 9y i Wo p ’ 1 1 y a y y i 3y 


Thus equation (84) becomes 


w 


U*A [1 

>1 (J 3C. 

j- + (L - 1) X) h. i 

Vp-Ve^ 

y i 1 3y 


(89) 
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It is seen that, when the Lewis number is unity, the heat- transfer rate is proportional 
to the enthalpy gradient at the edge of the equilibrium layer. Rather than setting the 
Lewis number equal to 1 at this point, the term explicitly displaying the gradient of 
the mass fraction is retained. 

The heat-transfer equation may be transformed into the rj coordinate using 
equations (57), (64), and (65) 


-«w* 


'Ve u e r T 

ISF 2 ? 


EC. h. „Z r . 
i l, e l, E i, 


: E ' + < L - - i 1 


(90) 


where the Prandtl number P = S 2 — 2 — is also not set equal to unity at this point. 

E 

Using equation (77), g’ „ and Z’. _ are eliminated yielding 


E 


i,E 


"^w = *e u e C H ( J e " *e) 


i + G K p) /l - 1) ? h ‘’ E ( C ‘ ,e C ‘ ,E ) 
\ L / I e -I E 


(91) 


where 


Pa u r 
_ e^e e 

p e u e C H ~ 


(2§) v “G(°o, P) 


(92) 


It is evident that, if the Lewis number is set equal to 1, the term in brackets 
in equation (91) becomes unity. For the case of L = 1, two implicit effects of mass 
injection on the heat flux are evident. First, mass injection causes an increase in 
G(°°,l) and a resultant decrease in the Stanton number C„. Secondly, variations in 

the enthalpy of the equilibrium layer, I E , will be caused by species variations. It is 

seen also that, when L > 1, an additional term is obtained which increases the heat 
transfer due to increased diffusion of chemical potential energy. 


Heat Balance on Fluid- Solid Interface 

In order to determine the net heat transferred from the boundary layer to the 
surface, it is necessary to examine the various energy fluxes crossing the gas-solid 
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interface. The following sketch schematically indicates the components of the energy 
balance at the interface. 



It is assumed that the gases flowing over the surface are nonradiating and trans- 
parent to the surface radiation q ,. The heat balance may be written 


“s * ■ ‘■rad + ril N h N 2 , E + rtl C h C, E (s) ‘ ,f, E h E < 8 ) < 93) 


The enthalpy of carbon in the solid phase is related to the enthalpy of carbon in 
the gaseous phase at the same temperature by the heat of sublimation of carbon as 
follows 


h c ( s ) + L vc = h c (g) (94) 


Therefore, combining equations (93) and (94), the heat flux to the surface becomes 


*■» * -*w - ‘‘rad + ri, N h N 2 , E + Iil C h C, E fe > ' *C L vc ' , *E h E (g) < 95 > 

The mass flux nig can be eliminated from equation (95) by using equation (10) to 
obtain 

- -«w - Vad + I *N h N 2 , E + “c h c, E (s) ' ™C L vc ' (™ 1 N + *c) h E (g) (96) 
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Finally, by using equation (91), the surface heat flux can be expressed as 


% 


" P e U e C H^e ' l E + G(~ S) ( L L *) £ h i, E ( C i, e " C i, e) " P 0 u *C R 
m / . m , \ m p 

* WiT i\- E ‘ hE<g) ) + Wii ( hc ' e(s) ' “ E<g) ) ' vS 


VC 


(97) 


It is seen from the enthalpy terms that the injection of nitrogen, in addition to 
those effects discussed previously, modifies the net energy flux to the surface. 


Stagnation- Point Problem 

Up to this point the analysis has been applicable to flat plate and stagnation-point 
boundary -layer problems, and to other configurations in the context of local similarity 
and negligible pressure gradient. Inasmuch as the stagnation -point problem is of 
greatest interest in current technology, this case is developed in more detail here. 

The velocity in the vicinity of a stagnation point is of the form 



(98) 


where 



is evaluated using the familiar Newtonian velocity gradient formula 



(99) 


From equation (58), 


s can be evaluated as 


s 



g2(j/+l) 

2(r + 1) 


( 100 ) 
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Equation (92) then yields 


P u C TT = 
e e H 



(v 


1/2 


+ 1 ) 
GK P) 


(101) 


The mass fluxes of nitrogen and carbon across the gas -solid interface are then 


. / *N \/„ \l/2 ( du e\ (v + l) l/2 

mN = y e u e C Hy ( Ve ) ' dS 4p G <°°’ P > 

A - /b An \ / p u \l/2(^e\ l/2 (^ + 1) 1/2 
C \ fi e U e C Hj P e ) ' dS 4p G <~ P > 


( 102 ) 


(103) 


The convective and diffusive heat fluxes at the wall for Lewis and Prandtl num- 
bers of 1 are given by 


. _ / \l/ 2 / du e\ ( y + l) 1 ^ 2 / _ \ 

q w “ ( P e P e) \dS/ G(°°, 1) \e ~ *Ey 


(104) 


and the heat flux penetrating the surface of the material is 


-/ 0 /du\^ 2 / i\i/2 r 

rf = (n n \^-/ ^ ( ®) ( v + 1) T _ 

q S ( P e p e) \dS/ gp GK 1) |^e 


m 


N 


Ie + \ Ve C H y * N 2’ e 


m 


+ B - 


N 


P e u e C H/ C ’ E 


V w fe> ' B 


m 


N 


L - Bh„(g) - 


rad 


p e u e G H/ VC E '*' p e u e C H 


(105) 


The mass flux of carbon across the interface is proportional to the rate of sur- 
face recession since 


• „ dT 

“t " P C dt 


(106) 
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where r is the thickness of the carbon material at the stagnation point. The surface 
ablation rate is then determined from 


dr _ _1_ 
dt = p c 



(, + D l/2 

G(°°, 1) 


(107) 


The thickness of carbon ablated during an arbitrary heating period is given by 
integration of equation (107) 


T = 



m 


B 


N 


P u C TT 
'e e H 



fr + 1) 
“gFT 


1/2 


1 ) 


dt 


(108) 


RESULTS AND DISCUSSION 


Calculations have been performed to demonstrate the accuracy of the approxi- 
mation of the present work in comparison with existing results for the case of no ni- 
trogen injection and to evaluate the effects of nitrogen injection. All calculations 
pertain to the stagnation-point boundary layer. The method of determination of 
boundary-layer-edge properties and composition has been outlined. Properties and 
species mass fractions, so obtained, are tabulated in table I. It is seen that at these 
flight conditions oxygen is completely dissociated and significant nitrogen dissociation 
has occurred. 

TABLE I. - FREE-STREAM AND STAGNATION- POINT CONDITIONS 

[Velocity = 19 700 ft/ sec] 


■ 1 

Case 

Alti- 

tude, 

kft 

V 

Btu/lb 

Pe’ 

atm 

T , 
e’ 

°K 

M /M 
e/ a 

P e' 

min 3 

C N, e 

C N 2 ,e 

C 0, e 

C °2’ e 

a 

115 

7864 

2.985 

6750 

1.379 

7.070 x 10" 3 

0. 19 

0. 58 

0.23 

<10" 3 

b 

150 

7875 

. 605 

6220 

1.402 

1. 530 

.21 

. 56 

.23 

<10" 3 

c 

180 

7877 

. 2026 

5900 

1.418 

. 535 

.21 

. 56 

.23 

<10“ 3 

d 

200 

7872 

.0966 

5750 

1.425 

.260 

.24 

| 

. 53 

. 23 

<10" 3 

e 

220 

7862 

.0458 

5560 

1.431 

. 126 

.25 

. 52 

.23 

<10" 3 
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Table II presents calculations of enthalpy, mass injection parameters, and 
species mass fractions in the equilibrium layer for a wall temperature of 2000° K. 
Equilibrium constants for equations (2) to (4) and species properties in the equilibrium 
layer were obtained from reference 13. In the diffusion -controlled regime it is found 
that only molecular nitrogen and carbon monoxide exist in significant amounts in the 
equilibrium layer. This is in agreement with the findings of reference 3 for no nitro- 
gen injection. At 2000° K, for instance, the ratios of the partial pressure of other 

species to the total pressure were found to be on the order of 10 ^ or less. It was 
also found from equation (43) that the relative proportions of nitrogen and carbon mon- 


m. 


oxide in the equilibrium layer depend only on the nitrogen injection parameter — ■ j-L- 

• p u C T 

m_ ^ ^ i 


Increasing values of — , of course, reduce the mass fraction of carbon monoxide 

p e e H 


U V TT 

e e H 


at the wall. 


TABLE H. - EQUILIBRIUM LAYER PROPERTIES 

[T w = 2000° kj 


“n 

C CO,E 

Cn 2 ,e 

m c 

B 

J E’ 

Btu/lb 


p e u e C H 

0 

0. 344 

0. 656 

0.173 

0. 173 

284 

.5 

. 242 

. 758 

.173 

. 673 

455 

1.0 

. 186 

. 814 

. 173 

1. 173 

548 

2.0 

. 128 

. 872 

. 173 

2. 173 

646 

3.0 

.097 

.903 

. 173 

3. 173 

698 ! 

4.0 

.078 

.922 

. 173 

4.173 

731 

6.0 

.056 

.944 

. 173 

6. 173 

769 


It can be demonstrated that the carbon mass flux parameter 
pendent of the nitrogen mass injection parameter as shown in table 


is inde- 


m c 

p e U e C H 
n. Elimination of 
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B + 1 in equation (23) using equation (25) allows expression of the carbon mass flux 
parameter in terms of boundary -layer -edge oxygen mass fraction only 


m c 3 _ 

p e u e C H = ?C °’ e 


(109) 


m, 


While 


p e u e C H 


is invariant with in- 


creases in the nitrogen injection parameter, 
it is evident that the carbon mass loss is 
reduced due to the decreased Stanton num- 
ber resulting from the increased total mass 
injection. The reduction in carbon mass 

rfu 


loss with increasing 


P u 
'e e 


N 

u 


and increas- 


H 


ing altitude is demonstrated in figure 2. 
Increasing the altitude reduces the carbon 
mass loss due to the decreasing flux of 
oxygen available from the free stream for 
reaction. 



Figure 2. - Effects of nitrogen injection 
and altitude on carbon mass loss. 


As discussed previously, perhaps the 
most accurate and sophisticated analysis for 
the ablation of graphite appears in refer - 

ence 3. For the special case of n ,, n = 0, 

P e u e^H 

the results of the present analysis are com- 
pared with the carbon mass loss obtained 
from reference 3. It is seen that the pres- 
ent results are consistently lower by about 
10 to 15 percent. Differences of this mag- 
nitude are not surprising in view of the vast 
differences between the present model and 
that of reference 3. 


The variation of total mass flux from 
the equilibrium layer is shown in figure 3 
as a function of the nitrogen injection pa- 
rameter. It is seen that the total mass flux 
increases rapidly for small values of the 
injection parameter. At larger values of 


m 


N 


p e u e C H 


the increase is less marked. This 


behavior results from the fact that C 


H 



0 1 2 3 4 5 


P«u* C H 

Figure 3. - Total mass loss rate. 
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decreases rapidly for large values, of mass injection, resulting in small increases in 


total mass flux as 


m 


Pe u e 


N 


H 


is increased. 



Altitude, ft 

Figure 4. - Aerodynamic heating rate. 


A comparison of the aerodynamic 
heat flux with carbon ablation from the 
present analysis and reference 3 is shown 
in figure 4. Differences of about 10 per- 
cent are evident. In view of the close cor- 
respondence of the present results to those 
of the sophisticated work of reference 3, it 
is felt that considerable confidence may be 
placed in the present analysis. It is also 
evident from figure 4 that nitrogen injection 
provides substantial reductions in heat flux. 
This is attributable primarily to large 
reductions in the Stanton number due to 
blowing. 

The heat flux absorbed by the graph- 
ite surface as calculated by equation (97) 
is shown in figure 5 as a function of the 
injection parameter for an altitude of 
115 000 feet. It is seen that nitrogen in- 
jection strongly reduces the heat flux 
absorbed. 



0 12 3 4 5 6 

<>« “e C H 


Figure 5. - Surface heat flux absorbed. 


The results presented have all as- 
sumed that an unidentified mechanism con- 
trolled the surface temperature at a fixed 
value. Usually, at a given flight condition, 
the wall temperature will be determined by 
the rate of radiation cooling, mass loss, 
and internal conduction. It is of interest to 
estimate the nitrogen injection required to 
maintain the wall at a specified wall tem- 
perature. In the case of steady ablation the 
surface temperature is determined by a 
balance of the net heat flux crossing the 
solid-gas interface and of the enthalpy trans- 
port of carbon and nitrogen from the initial 
temperature T in 


«h 



'p, C 


(s) dT + m 




( 110 ) 
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This expression is plotted in figure 5 for T in = 300° K and T = 2000° K. The 
intersection of the two curves shows that a wall temperature of 2000° K is obtained at 

the listed flight conditions with a nitrogen mass injection parameter = 3. 5. 

p e u e C H 

Figures 2 and 3 show that this corresponds to graphite and nitrogen expenditures of 

o 

0. 0029 and 0. 056 lb/ft -sec, respectively. 


A widely used parameter in evaluating the merit of a heat -protection system is 
the effective heat of ablation. It is a measure of the thermal energy accommodated per 
unit mass loss of the heat -protection system. The effective heat of ablation as con- 
ventionally defined is 


h = 
eff 



( 111 ) 


where 



B. 


is the aerodynamic heat transfer without mass injection. 


Using 


equations (103) and (104) with I- = I., the effective heat of ablation becomes 



The effective heat of ablation is plot- 
ted in figure 6 as a function of the enthalpy 

difference I -I A , where I. is the en- 
e A A 

thalpy of air at wall temperature and pres- 
sure. Again comparison with the results of 
reference 3 shows satisfactory agreement 
for the case of no nitrogen injection. It is 
seen that, as judged by h^, substantial 

gains are made in ablation efficiency as a 
result of mass injection. This result dem- 
onstrates that, if mass loss of graphite is 
critical, nitrogen injection may be employed 
to reduce surface ablation substantially. 
Clearly, other parameters, which show less 
favorable comparisons, may be defined to 
represent the ablation efficiency. Care 



( 112 ) 



o 2468 10X10 3 . 

Ie-I A ,Btu/lb 


Figure 6. - Effective heat of ablation. 
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should be taken, therefore, in judging heat-protection-system efficiency based on 
parameters such as h^. 


CONCLUDING REMARKS 


The effects of nitrogen injection on the heat and mass transfer of an ablating 
surface have been analyzed based on a model which considers all chemical reactions to 
occur in an infinitesimal region adjacent to the surface. Diffusion of oxygen through a 
frozen boundary layer and chemical equilibrium in the equilibrium layer have been as- 
sumed to control the rate of mass loss of carbon. Equations have been derived for the 
species mass fractions and mass injection parameters as well as for the heat and mass 
transfer rates. It has been demonstrated that: 

(1) The carbon mass injection parameter is independent of nitrogen injection in 
the diffusion-controlled regime and, in fact, depends only on the mass fraction of oxy- 
gen at the boundary -layer edge. 

(2) The only species present at the wall in significant quantities are nitrogen and 
carbon monoxide. Injection of nitrogen increases the wall mass fraction of nitrogen at 
the expense of the carbon monoxide mass fraction. 

(3) The injection of nitrogen affects heat transfer by reducing the heat transfer 
coefficient as a result of increased blowing into the boundary layer and through modi- 
fications of the enthalpy of gases at the wall. 

(4) When minimization of carbon surface recession is important, injection of 
nitrogen can substantially reduce the carbon mass loss. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, March 23, 1967 
914-50-10-01-72 
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